clear all
adjlink = csvread('D:\shale_job\neighborCounty.csv');
adjmax = size(adjlink, 1);
countyNo = 254;
adjmatrix = zeros(countyNo);
for i = 2:adjmax
    for j = 1:countyNo
        if adjlink(i,j) ~= 0
            ind = adjlink(i,j);
            adjmatrix(ind,j) = 1;
        end
    end
end
adjmax_rstd = adjmatrix./repmat(sum(adjmatrix,2),1,254);
    